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^ \ There are good reasons to consider models of low-energy supersymmetry 

^ \ with very light photinos and gluinos. In a wide class of models the lightest R- 

^ ■ n 

O ' odd, color-singlet state containing a gluino, the i? , has a mass in the 1-2 GeV 

in : 

^ ' range and the slightly lighter photino, 7, would survive as the relic i?-odd 

Q-i| species. For the light photino masses considered here, previous calculations 

o ■ 

?-j \ resulted in an unacceptable photino relic abundance. But we point out that 

^ ' processes other than photino self-annihilation determine the relic abundance 

> ' 

I when the photino and R^ are close in mass. Including R^ < — > 7 processes, 

I we find that the photino relic abundance is most sensitive to the R^-to-'y 

mass ratio, and within model uncertainties, a critical density in photinos 
may be obtained for an R^-to-'j mass ratio in the range 1.2 to 2.2. We 
propose photinos in the mass range of 500 MeV to 1.6 GeV as a dark matter 
candidate, and discuss a strategy to test the hypothesis. 
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I. INTRODUCTION 



In this paper we study the early-Universe evolution and freeze out of light, long-lived 
or stable, i?-odd states, the photino 7 and the gluino g.^ In the type of models we 
consider, the photino should be the relic -R-odd particle, even though it may be more 
massive than the gluino. This is because below the confinement transition the gluino 
is bound into a color-singlet hadron, the whose mass (which is in the 1 to 2 GeV 
range when the gluino is very light [jl],^) is greater than that of the photino. Including 
previously neglected reactions associated with the gluino (more precisely, associated with 
the we find that light photinos may be cosmologically acceptable; indeed they are 
an attractive dark-matter candidate. 

In the minimal SUSY model, the mass matrix of the charged and neutral SUSY fermions 

(gauginos and Higgsinos) are determined by Lagrangian terms involving the Higgs chiral 

superfields. Hi and H2, and the SU(2) and U(l) gauge superfields, and B, plus soft 

supersymmetry breaking terms. This leads to a neutralino mass matrix in the basis 

{B, W^, Hf, H^) of the form 

/ Ml —mz cos P sin 9w f^z sin P sin 9]y \ 

M2 mz cos /? cos 6w — "^z sin [3 cos 6y/ 

—mz cos (3 sin 9w mz cos (3 cos 9w — /i 

\ sin/5 sin — sin/5cos6'vi/ — /x / 

Here mz is the mass of the Z-boson, 9w is the Weinberg angle, /x is the coefficient 

of a supersymmetric mixing term between the Higgs superfields, and tan (3 is the ratio 

of the vacuum expectation values of the two Higgs fields responsible for electroweak 

symmetry breaking. The suSY-breaking masses Mi and M2 are commonly assumed to 

^i?-parity is a multiplicative quantum number, exactly conserved in most SUSY models, under which 
ordinary particles have i? = +1 while new "superpartners" have i? = — 1. Throughout this paper we 
will assume that i?-parity is exact so the lightest i?-odd particle is stable. 
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be of order mz or larger, and if the SUSY model is embedded in a grand unified theory, 
then 3M1/M2 = 5ai/a2. 

The terms in the Lagrangian proportional to Mi and M2 arise from dimension-3 
SUSY-breaking operators. However such SUSY-breaking terms are not without problems. 
It appears difficult to break SUSY dynamically in a way that produces dimension-3 terms 
while avoiding problems associated with the addition of gauge-singlet superfields In 
models where SUSY is broken dynamically or spontaneously in the hidden sector and there 
are no gauge singlets, all dimension-3 SUSY-breaking operators in the effective low-energy 
theory are suppressed by a factor of ($)/mp;, where ($) is the vacuum expectation value 
of some hidden-sector field. Thus, dimension-3 terms effectively do not contribute to the 
neutralino mass matrix. This would imply that at the tree level the gluino is massless, 
and the neutralino mass matrix is given by Eq. (|l]) with vanishing (00) and (11) entries. 
However, a non-zero gluino mass, as well as non-zero entries in the neutralino mass 
matrix are generated through radiative corrections such as the top-stop loop, and for the 
neutralinos, "electroweak" loops involving higgsinos and/or winos and binos. 

The generation of radiative gaugino masses in the absence of dimension-3 SUSY break- 
ing was studied by Farrar and Masiero [^.0 From Figs. 4 and 5 of that paper one sees 
that as Mq, the typical SUSY-breaking scalar mass, varies between 100 and 400 GeV, the 
gluino mass ranges from about 700 to about 100 MeV,0 while the photinoQ mass ranges 

^See also 1^ for general formulae. Earlier studies ^,0] of radiative corrections when tree level gaugino 
masses are absent included another dimension-3 operator, the so-called "A term," and did not consider 
the electroweak loop contributions to the neutralino mass matrix. They also assumed model-dependent 
relations between parameters. 

•^Actually, larger values of Mq are not considered in order to keep the gluino mass greater than about 
100 MeV. Otherwise an unacceptably light pseudoscalar meson would be produced [Q. 

*Upon diagonalization of the mass matrix, the physical neutralino states are a linear combination 
of B'^ , W^, H'^, and iJj- When the gaugino submatrix elements are small, the lightest neutralino is 
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from around 400 to 900 MeV, for /i > 40 GeV. This estimate for the photino mass should 
be considered as merely indicative of its possible value, since an approximation for the 
electroweak loop used in Ref. is strictly valid only when fi or Mq are much larger than 
m^y. The other neutralinos are much heavier, and the production rates of the photino 
and the next-lightest neutralino in Z° decay are consistent with LEP bounds @]. 

Using the results of Ref. @| , but additionally restricting parameters so that the cor- 
rect electroweak symmetry breaking is obtained, Farrar 0] found Mq ~ 150 GeV and 
estimated the lifetime. This allowed completion of the study of the main phenomeno- 
logical features of this scenario, which was begun in Ref. 0]. The conclusion is that 
light gluinos and photinos are quite consistent with present experiments, and result in a 
number of striking predictions . However models with light gauginos have been widely 
thought to be disallowed because it has been believed that the relic density of the lightest 
neutralino, usually referred to as the LSP,0 exceeds cosmological bounds unless i?-parity 
is violated [§-0 . 

In this paper we point out that previous considerations of the relic abundance have 
neglected the rather important interplay between the photino and the gluino which can 
determine the final neutralino abundance if the photino and gluino are both light, as 
they must be in models without dimension-3 explicit suSY-breaking terms. We find 
that when gluino-photino interactions are included, rather than being a cosmological 
embarrassment, these very light photinos are an excellent dark matter candidate. In this 
paper we discuss the decoupling and relic abundance of light photinos, and the sensitivity 
of the result upon the parameters of the SUSY models. 

a linear combination of and i?° that is almost identical to the SU(2) x U(l) composition of the 
photon, and thus is correctly called "photino." 

^In this scenario, lsp is an ambiguous term: the gluino is lighter than the photino, although the 
photino is lighter than the . A more relevant term would be lrocs — lightest R-odd color singlet. 
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For the light masses studied here, freeze-out occurs well after the confinement tran- 
sition so the physical states must be color singlets. Since g is not a color singlet, below 
the confinement transition the relevant state to consider is the lightest color-singlet state 
containing a gluino, which is believed to be a gluon-gluino bound state known as the R^. 
The other light -R-odd states are more massive than these, and decay to the two light 
ones with lifetimes much faster than the expansion rate at freeze out. The only other 
possible state of interest is the 5", which is the lightest i?-odd baryon, consisting of the 
color-singlet, fiavor-singlet state udsg [jl],|^,[lO| . The masses for the states we consider will 
be assumed to be in the range [0,^ 

^(gluino) : m~ = 100 - 600 MeV 

7(photino) : m = 100 MeV - 1400 MeV 

R%gg) : M = 1 - 2 GeV 

S^udsg) : M50 = 1.5-2 GeV. (2) 

Since it is the lightest color-singlet i?-odd state, the 7 is stable, and decays to a final 
state consisting of a photino and typically one meson: R^ — > 77?; 777, etc. The lifetime 
is very uncertain, but probably lies in the range 10~^ to 10~^°s, or even longer p|. 

The reaction rates that determine freeze out will depend upon the 7 and i?" masses, 
the cross sections involving the 7 and and possibly the decay width of the R^ as well. 
In turn the cross sections and decay width also depend on the masses of the 7, g and 
as well as the masses of the squarks and sleptons. We will denote the squark/slepton 
masses by a common mass scale (expected to be of order 100 GeV). Even if the 
masses were known and the short-distance physics specified, calculation of the width 
and some of the cross sections would be no easy task, because one is dealing with light 
hadrons. Fortunately, our conclusions are reasonably insensitive to individual masses, 
lifetimes, and cross sections, but depend crucially upon the R'^-to-j mass ratio. 
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When we do need an explicit value of the photino mass m or the masses of squarks 

and sleptons, we will parameterize them by the dimensionless ratios 

_ m _ 

= 800 MeV' = 100 GeV " 

Although there are several undetermined parameters in our calculation, as mentioned 
above, the most important parameter will be the the ratio of the mass to the 7 mass: 

r = —. 4 

m 

This is by far the most crucial parameter, with the relic abundance having an exponen- 
tial dependence upon r. We find that limits to the magnitude of the contribution to 
the present mass density from relic photinos requires r < 2.2,0 while r must be larger 
than about 1.2 if the photino relic density is to be significant. This narrow band of r 
encompasses the large uncertainties in lifetimes and cross sections. If the mass ratio is 
between about 1.6 and 2, then light-mass relic photinos could dominate the Universe 
and provide the dark matter with fi- ~ 1. 

In the concluding section we explore the proposal that light photinos are the dark 
matter, and discuss possibilities for testing the idea. We lay the groundwork for this 
suggestion in the next section as we develop a new scenario for decoupling and freeze- 
out for the photinos and gluinos. In Section III we consider the cross sections and 
lifetimes used in Section IV to calculate the reaction rates relevant for the determination 
of the freeze-out abundance of the photinos (and hence Q-h"^). In Section V we compare 
the reaction rates to the expansion rate and estimate when photinos decouple. 



^Or else i?-parity must be violated so the photinos decay. 
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II. SCENARIO FOR PHOTINO/GLUINO FREEZE OUT 



The standard procedure for the calculation of the present number density of a thermal 
relic of the early-Universe is to assume that the particle species was once in thermal 
equilibrium until at some point the rates for self-annihilation and pair-creation processes 
became much smaller than the expansion rate, and the particle species effectively froze 
out of equilibrium. After freeze out, its number density decreased only because of the 
dilution due to the expansion of the Universe. (For a discussion, see Ref. [0.) 

Since subsequent to freeze out the number of particles in a co-moving volume is 
constant, it is convenient to express the number density of the particle species in terms 
of the entropy density, since the entropy in a co-moving volume is also constant for most 
of the history of the Universe. The number density-to-entropy ratio is usually denoted 
by y. If a species of mass m is in equilibrium and non-relativistic, Y is simply given in 
terms of the mass-to-temperature ratio x = m/T as 



where g is the number of spin degrees of freedom, and (yf* is the total number of relativistic 
degrees of freedom in the Universe at temperature T = m/x. Well after freeze-out Y{x) 
is constant, and we will denote this asymptotic value of Y as 1^00- 

If self annihilation determines the final abundance of a species, Yoo can be found by 
integrating the Boltzmann equation (dot denotes d/dt) 



where n is the actual number density, n^Q is the equilibrium density, H is the expansion 
rate of the Universe, and (|f Icr^i) is the thermal average [0,0 of the annihilation rate. 

There are no general closed-form solutions to the Boltzmann equation, but there 
are reliable, well tested approximations for the late-time solution, i.e., Fqo- Then with 



Yeq{x) = 0.145((7/(?,)x3/2exp(-x) 



(5) 




(6) 
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knowledge of Yoo, the contribution to flh'^ from the species can easily be found. Let us 
specialize to the survival of photinos assuming self-annihilation determines freeze out. 

Calculation of the relic abundance involves first calculating the value of x, known 
as Xf, where the abundance starts to depart from the equilibrium abundance. Using 
standard approximate solutions to the Boltzmann equation [0 givesQ 

Xf = ln(0.0481mpim(Jo) - 1.5 ln[ln(0.0481mpimao)] , (7) 

where we have used g = 2 and g^, = 10, and parameterized the non-relativistic annihila- 
tion cross section as (|f |(Ta) = ctqx^^. In anticipation of the results of the next section, 
we use (To = 2 X lO^^^filfig'^mb, and we find Xf ~ 12.3 + ln(/ig//i|). The value of Xf 
determines Foo^ 

= ^ ^ 7.4 X 10-^8 4. 8 

mpimao 

Once is known, the present photino energy density can be easily calculated: 
p;^ = mn~ = 0.8 /isGeV ■ Foo 2970 cm"^. When this result is divided by the critical 
density, pc = 1.054/1^ x 10^^ GeVcm^^, the fraction of the critical density contributed 
by the photino is fi-Zi^ = 2.25 x lOVs^oo- For Y^o in Eq. (|), fi-Zi^ = 167p^V|- 

The age of the Universe restricts ^l-h"^ to be less than one, so for = 1, the photino 
must be more massive than about 10 GeV if its relic abundance is determined by self- 
annihilation. 

But in this paper we point out that for models in which both the photino and the 
gluino are light, freeze-out is not determined by photino self annihilation, but by 7-/2° 
interconversion. The basic point is that since the has strong interactions, it will 
stay in equilibrium longer than the photino, even though it is more massive. As long 



^Freeze-out aficionados will notice that we use the formulae appropriate for p-wave annihilation 
because Fermi statistics requires the initial identical Majorana fermions to be in an L = 1 state 
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as 7 ^ interconversion occurs at a rate larger than H, then through its interactions 
with the the photino will be able to maintain its equilibrium abundance even after 
self-annihilation has frozen out.0 

Before we demonstrate that this scenario naturally occurs for the types of photino 
and masses expected, we must determine the cross sections and decay width of the 
reactions involving the photino and the gluino. 



III. CROSS SECTIONS AND DECAY WIDTH 



In this section we characterize the cross sections and decay width required for the 

determination of the relic photino abundance, and also discuss the uncertainties. We 

should emphasize that all cross sections are calculated in the non-relativistic (N.R.) 

limit, and by (■ ■ •) we imply that the quantity is to be evaluated as a thermal average 

|]T^,|TB|. In the N.R. limit a temperature dependence to the cross sections enters if the 

annihilation proceeds through a p-wave, as required if the initial state consists of identical 

fermions . For p-wave annihilation, at low energy the cross section is proportional to 

f^, where |f | is the relative velocity of the initial particles. The thermal average reduces 

^Actually, interconversion can also play an important role in determining the relic abundance of 
heavier photinos. When the photino is more massive and freeze-out occurs above the confinement phase 
transition, the analysis is similar to the one here; in fact simpler because perturbation theory can be 
used to compute the relevant rates involving gluinos and photinos. Since the qualitative relation between 
interconversion and self-annihilation rates is independent of whether the gluino is free or confined in an 
i?'^, one can get a crude idea of the required gluino-photino mass ratio, r, just by using the analysis in 
this paper and scaling the results to the value of /ig of interest. We concentrate on the light gaugino 
scenario because it is attractive in its own right, and also because it naturally produces r in the right 
ballpark[^. In a conventional SUSY-breaking scheme fine-tuning is generally necessary to give r the right 
value for the interconversion mechanism to play an important role. 
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to replacing by QT/m, where m is the mass of the particle in the initial state. 

We now consider in turn the cross sections and width for the individual reactions 
discussed in the previous section. 

A. Self-annihilations and co-annihilation 

The first type of reactions we will consider are those which change the number of 
i?-odd particles. 

BPR^ — » X : We will refer to this process as self-annihilation. At the constituent 
level the relevant reactions are g+g —>■ g+g and g+g —>■ q+q, which are unsuppressed by 
any powers of Mj, and should be typical of strong interaction cross sections. In the N.R. 
limit, we expect the R^R^ annihilation cross section to be comparable to the pp cross 
section, but with an extra factor of f ^, accounting for the fact that there are identical 
fermions in the initial state, so annihilation must proceed through a p-wave.^ There is 
some energy dependence to the pp cross section, but it is sufficient to consider {\v\ajiOjio) 
to be a constant, approximately given by 

{\vWroro) ^ lOOf^mb = QOOx-^r-^mh, (9) 

where we have used for the relative velocity = QT/M = 6/{rx), with x = m/T. 

We should note that the thermal average of the cross section might be even larger 
if there are resonances near threshold. In any case, this cross section should be much 
larger than any cross section involving the photino, and will ensure that the R^ remains 
in equilibrium longer than the 7, greatly simplifying our considerations. 

77 ^ X: In photino self-annihilation at low energies the final state X is a lepton- 

^In general the result is not so simple. For instance, in addition to the term proportional to , 
the cross section also involves a term proportional to the square of the masses of the initial and final 
particles. 



9 



antilepton pair, or a quark-antiquark pair which appears as hght mesons. The process 
involves the t-channel exchange of a virtual squark or slepton between the photinos, 
producing the final-state fermion-antifermion pair. In the low-energy limit the mass 
Mgr of the squark/slepton is much greater than ^/i, and the photino-photino-fermion- 
antifermion operator appears in the low-energy theory with a coefficient proportional to 
ef/Mj^, with Cj the charge of the final-state fermionj^ Also, as there are two identical 
fermions in the initial state, the annihilation proceeds as a p-wave, which introduces 
a factor of v'^ in the low-energy cross section [Q. The resultant low-energy photino 
self-annihilation cross section is 0,^,14,1^: 



2 2 

m V 



(|t;|a:^) = 87ra|,,Egf^y ^2.0x10-11 a;-i mb, (10) 

where we have used for the relative velocity v"^ = 6/x with x = m/T, and is the 
magnitude of the charge of a final-state fermion in units of the electron charge. For the 
light photinos we consider, summing over e, fi, and three colors of u, d, and s quarks 
leads to J2i if = 8/3. 

7/?" — > X: This is an example of a phenomenon known as co-annihilation, whereby 
the particle of interest (in our case the photino) disappears by annihilating with another 
particle (here, the i?°). Of course co-annihilation also leads to a net decrease in i?-odd 
particles. 

In all processes involving the photino--R° interaction, the leading tree-level short- 
distance operator containing g and 7 is A~A~g|gj + h.c, with coefficient eqiQs/M^'^. For 
three light quarks, J2i Qi = 2. Thus we can estimate the cross section for jR^ — > X in 
terms of the 7 self annihilation cross section: 

, , ^5 4 2 M 3 , , , , 

ilvW^Ro) ^ 2 ' 11 

CiEM 3 8/3 m f ^ " 



^°The electric charge e and the strong charge gs are to be evaluated at a scale of order M~, so in 
numerical estimates we use asm — 1/128 and as — 0.117. 
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where the ratio of a's arises because the short-distance operator for co-annihilation is 
proportional to efgg rather than ef, the second factor is the color factor coming from 
the gluino coupling, and the third factor comes from the ratio of J2i Qi / J2i Qi for the 
participating fermions. We have replaced appearing in Eq. (|IU]) by mM, although 
the actual dependence on m and M may be more complicated. Finally, the annihilation 
is s-wave so there is no suppression as in photino self-annihilation. 

Although the short- distance physics is perturbative, the initial gluino appears in a 
light hadron, and there are complications in the momentum fraction of the carried 
by the gluino and other non-perturbative effects. For our purposes it will be sufficient 
to account for the uncertainty by including in the cross section an unknown coefficient 
A, leading to a final expression 

mb. (12) 

It is reassuring that if one estimates {\v\a~fio) starting from {\v\aROi^o) a similar result is 
obtained. We find that co-annihilation will not be important unless A is larger than 10^ 
or so, which we believe is unlikely. 

B. 'j—R^ inter conversion 

In what we call interconversion processes, there is an i?-odd particle in the initial as 
well as the final state. Although the reactions do not of themselves change the number 
of i?-odd particles, they keep the photinos in equilibrium with the R^s, which in turn 
are kept in equilibrium through their self annihilations. 

i?" —>■ 77r: R^ decay can occur via, e.g., the gluino inside the initial R^ turning into an 
antiquark and a virtual squark, followed by squark decay into a photino and a quark. In 
the low-energy limit the quark-antiquark-gluino-photino vertex can be described by the 
same type of four-Fermi interaction as in co-annihilation. One expects on dimensional 
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V a; 



) ~ 1.5 X 10"^° r 



^l^s^A 



grounds a decay width Fq oc aEMdsM^ / The hfetime of a free gluino to decay to 



a photino and massless quark-antiquark pair was computed in Ref. |T^. However this 



does not provide a very useful estimate when the gluino mass is less than the photino 
mass. 

The lifetime for decay was studied in Ref. In an attempt to account for 
the effects of gluino-gluon interactions in the necessary for even a crude estimate 
of the lifetime, the following picture was developed: The is viewed as a state 
with a massless gluon carrying momentum fraction x, and a gluino carrying momentum 
fraction (1 — x)J^ having therefore an effective mass M\/l — x. The gluon structure 
function F{x) gives the probability in an interval x to x + dx oi finding a gluon, and 
the corresponding effective mass for the gluino. One then obtains the R^ decay width 
(neglecting the mass of final state hadrons): 



where Tq{M^ 0) is the rate for a gluino of mass M to decay to a massless photino, and 



f[y) = [(1 - y^){l + 2y- 7y^ + 20y^ - 7y^ + 2y^ + y^) + 24?/3(l - y + y'^) logy] contains 



the phase space suppression which is important when the photino becomes massive in 

comparison to the gluino. Modeling decay in a similar manner underestimates the 

lifetime by a factor of 2 to 4. This is in surprisingly good agreement; however caution 

should be exercised when extending the model to R^ decay, because kaon decay is much 

less sensitive to the phase-space suppression from the final state masses than the present 

case, since the range of interest will turn out to be r ~ 1.2 — 2.2. For r in this range, 

taking F{x) ~ 6x(l — x) following the discussion in Ref. [Q] leads to an approximate 

^^Of course there should be no confusion with the fact that in the discussion of the hfetime we use 
X to denote the gluon momentum fraction whereas throughout the rest of the paper x denotes ni/T . 
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(13) 
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behavior 

r^o_:^. = 2.0 X 10-^^ Hr) GeV [^4^^tB] , (14) 

where J^{r) = r^(l — r~^)^, and the factor B reflects the overall uncertainty. We believe 
a reasonable range for B is 1/30 < -B < 3. 

R^n < — > 77r: We will refer to these processes as photino-i?*^ conversion, since an 
initial (or 7) is converted to a final 7 (or The short-distance subprocess in this 
reaction is g + ^ — g + 7, again described by the same low-energy effective operator as in 
co-annihilation and decay. At the hadronic level the matrix element for R^tt —* 7X is 
the same as for R^j —>■ nX for any X, evaluated in different physical regions. Thus the 
difference between the various cross sections is just due to the difference in fluxes and 
final state phase space integrations, and variations of the matrix element with kinematic 
variables. Given the crude nature of the analysis here, and the great uncertainty in 
the overall magnitude of the cross sections, incorporating the constraints of crossing 
symmetry are not justified at present. 

We can point to one specific hadronic effect which we do not include but which is 
potentially important. It is likely that near threshold there is a resonance (the Rj^) which 
would increase the cross section by a factor of 4M|,/r^, where M/j is the mass and 
the width of the resonance. This complicates matters because neither the resonance's 
width nor its distance above threshold is known. If a resonance is important, it would 
also be necessary to perform the thermal average over the resonance in a more careful 



manner 



We will parameterize our uncertainty by including a factor C in the cross section to 
express the uncertainty due to hadronic physics and the possible existence of a resonance. 
Putting everything together, we obtain 

{\v\aRo^) ^ 1.5 X 10^^° r [f^l^s^C] mb. (15) 
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Table I: Cross sections and the decay width used in the calculation of the relic photino abundance. 
The dimensionless parameters fig and fis were defined in Eq. (^), and J-{r) was discussed below Eq. 
(p^). The coefficients A, B, and C reflect uncertainties involving the calculation of hadronic matrix 
elements. 



Process 


Cross section or width 


i?'' self annihilation: (|w|(T/jo/jo) 
7 self annihilation: 
co-annihilation: {\v\o'~j^o) 
R° decay: T^,,^-^ 
7 - i?° conversion: {\v\<jj^o^) 


600 x^^ mb 
2.0 X f 0"" x-^ [mIa^s''] ™b 
1.5 X 10-1° r [nlfig'^A] mb 
2.0 X 10-1"^ T{r) [mIms^-S] GeV 
1.5 X 10-1" r [filfig^C] mb 



We would expect C to fall in the range 1 < C < 10^. We will use detailed balance 
arguments which allow us to avoid using the inverse reaction, 77? R^tt. 

This completes the discussion of the lifetimes, cross sections, and their uncertainties. 
The results are summarized in Table I. 



IV. EARLY-UNIVERSE REACTION RATES 

To obtain an estimate of when the rates will drop below the expansion rate, we will 
assume all particles are in LTE (local thermodynamic equilibrium). In LTE a particle of 
mass m in the N.R. limit has a number density 



n 



-^-^(mTf expi-m/T) = -^^^(T/mf'm'expi-m/T). (16) 

Here g counts the number of spin degrees of freedom, and will be 2 for the and the 7. 

H (the expansion rate): Of course all rates are to be compared with the expansion 
rate. In the radiation-dominated Universe with g^: ~ 10 degrees of freedom 

H = l.mgl'^T^/mpi = 2.8 x lO-^^a;-^ Ifij] GeV. (17) 
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77 



X (photino self-annihilation): In the Boltzmann equation for the evolution of 



the 7 number density there are terms accounting for photino self-annihilation and photino 
pair-production from light particles in the plasma. Assuming the light annihilation 
products are in LTE, the terms are of the form 



n~ + 3Hn~ D 



(18) 



If we assume that the photino is in equilibrium, the self-annihilation and pair production 
terms are equal, and we may express the individual terms in the form 



n 



EQ 



(19) 



where the upper sign is for self-annihilation and the lower sign is for pair production. 



It is obvious that 



n^^ {\v\a~z) plays the role of a "rate" to be compared to H. If 



this rate is much greater than the self-annihilation/pair-production processes will 
ensure the photino is in equilibrium, while if the rate is much smaller than if, self- 
annihilation/pair-production cannot enforce equilibrium. 

Therefore we define an equilibrium photino annihilation rate by r(77 X) = 
rz^'^df |o"~). Using Eq. (|^) for the equilibrium abundance and the annihilation cross 
section discussed in the previous section, we find 

2 /T\^/2 , 2.0 X lO-i^mb 



r(77 X) 



GeV. 



(27r)3/2 

3.3 X 10"^2a;~^/2gxp(-x) 



X 



-1 



2 -4 



(20) 



R^R^ X {R^ self-annihilation): Determination of the equilibrium rate for i?'^-self- 
annihilation proceeds in a similar manner, yielding r{R^R^ ^ X) = n^? {\v\(Tj^oro): 



r(i?"/?" X) 



rp . 3/2 



M3exp(-M/T) 



240 X ^ r ^ mb 



(27r)3/2 \Mj " " ' ^ " " ' 0.4 mb GeV" 
99r^/2a;-^/2exp(-rx) [i^l] GeV. 



(21) 
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'J BP — > X (7 co-annihilation): In the Boltzmann equation for the evolution of the 7 
density will appear a term —njn)n~{\v\a~^o) ■ Therefore the equilibrium co-annihilation 
rate for the decrease of the 7 density is 

TijR'^X) = n^o^ (|^|a~^o) 

= 2.5 X 10-^^ r^/" exp(-rx) [jil/is^] , (22) 

where we have again assumed the particles in the process are in equilibrium. 

77r i?° (Inverse decay) : If the i?° decay products (in this case 7 and tt) are in 
equilibrium, then the Boltzmann equation for the evolution of contains a term 

fiRO + SHuro D -Tr3^~^ (uro - n^?) . (23) 

The first term on the rhs represents decay, while the second term represents "inverse 
decay." Since inverse decay turns a 7 into a there will be an "inverse decay" term 
in the equation for the evolution of the 7 number density: 

h~ + 3Hn~ D -r^o_~, n^R? . (24) 

The right hand side can be written as nS'^ [n^^ / n^'^) rR;)_^~^. Therefore the inverse 
decay rate in the evolution of the photino number density contributes a term 

rij^^R') = r^o^:^. (n^o^/r^^) = r^o^;^. (- j exp[-(M-m)/r] 

= 2.0 X 10"^^ J^(r) exp [-(r - l)x] [l^ll^s^B] GeV. (25) 

77r R^n (photino--R° conversion): It is easiest to obtain this term by first consid- 
ering the term in the equation for uro due to the reverse process and then using detailed 
balance: 

fiRp D -n^nRo{\v\aRO^) . (26) 
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Table II: The ratio of the equilibrium reaction rates to the expansion rate for the indicated reactions. 
Shown in [• • •] is the scaling of the rates with unknown parameters characterizing the cross sections and 
decay width. 



Process 


V/H 


7 self annihilation 




X 


1.2 X 10'^ a;-i/2 exp(-x) 






BP self annihilation 




X 


3.5 X 102°a;-V2 rV2 exp(-ra;) 




[ms] 


CO- annihilation 


7i?0- 


X 


8.9 X 10'^ a;^^ r^/^ exp(-ra:) 






inverse decay 


77r —^ 


R° 


7.1 X 10^ a;2r3/2jr(r) exp[-(r- 


l)x] 




7 - i?° conversion 


77r — > 


R°Tr 


9.6 X 10*5 a;V2r-5/2 exp[-(r - l)x 


exp(-0.175/z^^a;) 


r 3/2 -4/-,i 



Since the photino-i?° conversion process creates a 7 there is a similar term in n~ with the 
opposite sign. Now we can write this in a form to calculate the rate for 7 annihilation 

by 



h~ D -fiRO = n-;,njio{\v\aRO^} 



n~ 



(27) 



Assuming equilibrium as before, the rate keeping the 7 in equilibrium can be ex- 
pressed as 



r(77r ^ i?°7r) 



7 



= 2.7 X 10-^2 r^/^ x-^''^ exp(-0.175/X8 ^a;) 
X exp[-(r - 1)0;] \piT^itC] . 



(28) 



Of course it is the ratio of the reaction rates to the expansion rate that will be used 
to estimate photino freeze out. These ratios are given in Table II. 

There are two striking features apparent when comparing the magnitudes of the 
equilibrium reaction rates in Table II. The first feature is that the numerical factor in 
self annihilation is enormous in comparison to the other numerical factors. This simply 
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reflects the fact tliat i?" anniliilation proceeds tlirougli a strong process, wliile tlie otlier 
processes are all suppressed by a factor of 

The other important feature is the exponential factors of the rates. They will largely 
determine when the photino will decouple, so it is worthwhile to examine them in detail. 

The exponential factor in 7 self-annihilation is simply e~™/"^, which arises from the 
equilibrium abundance of the 7. It is simple to understand: the probability of one 7 
to find another 7 with which to annihilate is proportional to the photino density, which 
contains a factor of e"™/^ in the N.R. limit. 

The similar exponential factor in self-annihilation is also easy to understand. An 

must find another R'^ to annihilate, and that probability is proportional to e~*^/^ = 

Co-annihilation is also an exothermic process, so the only exponential suppression is 
the probability of a 7 locating the for co-annihilation, proportional to the equilibrium 
number density of R^, in turn proportional to e~^'^f^ = e~^^ 

In inverse-decay the exponential factor is e"'-'""^-'^ = ^-i^^-'^)/'^ , The number density 
of target pions is e"™'"/-^, so this factor is present. It is necessary for the 71-7 collision 
to have sufficient center-of-mass energy to create the This introduces an addition 
suppression of e-(^^-"^-™'r)/T_ Combining the two exponential factors gives the result in 
Table II. 

Finally, photino-gluino conversion involves two exponential suppression factors. The 
first, e"™'"/"^ = e~°-^''^'^8 ^ represents the suppression in the pion number density,^ and 
since the mass of the R^ is greater than the mass of the 7, there is an additional g-CA^-"^)/^ 
suppression. 

The factors of x and r originate from three places: a factor of comes from dividing 

^^At the temperatures of interest for decoupling, pions might be cheap, but they are not free. 
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X = m~/T 



r(7R0^X)/H 

r(77T^R0)/H 

r/H = 1 



r(77-X)/H 

r(77T^R07T)/H 

r(RORo^X)/H 



Fig. 1: Equilibrium reaction rates divided by H iov r = 1.25 and 1.5, assuming = /is = 1, and 
that the factors A = B = C = 1. The rates can be easily scaled for other choices of the parameters. 
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Fig. 2: The same as Fig. 1, but for r = 1.75 and 2. 
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the rates by H, factors of r and x arise for pre-exponential factors in the number density, 
and finally they may appear explicitly in the cross section or decay width. 

The equilibrium reaction rates divided by H are shown in Figs. 1 and 2 for r = 1.25, 
1.5, 1.75, and 2. In the figures we have assumed fig = fJ^s = ^ = B = C = 1. Using the 
information in Table II it is possible to scale the curves for other values of the parameters. 



V. ANALYSIS 

Rather than integrate a complete reaction network for the evolution and freeze out 
of the photinos, we will assume that the photinos remain in equilibrium so long as there 
is a reaction depleting the 7 abundance that is larger than H. We will then assume that 
as soon as the rate of the last such reaction drops below H, the photinos immediately 
freeze out, and the photino-to-entropy ratio is frozen at that value. We will call this 
approximation the "sudden" approximation. 

We can get some idea of the accuracy of the sudden approximation by considering a 
simple system involving only photino self-annihilation. As discussed in Section II, there 
is a well developed formalism for calculating the self-annihilation freeze-out of a N.R. 
species |lll. Using that formalism in Section II, Eq. (|), we found Fqo — 7.4 x 10"''. 



Now let us compute Yoo using the sudden approximation. From Fig. 1 or Fig. 2, we 
see that r(77 X) = H at x = 14.7, independent of r. We will denote by the value 
of X when T = H. Using the sudden approximation that the 7 is in LTE until x = x* and 
immediately freezes out would give a photino to entropy ratio of (again using 2 degrees 
of freedom and = 10) 

Yoo = YsQix,) = 0.145 (2/10)x^/2 exp(-x*) = 7 x 10"^ (using x, = 14.7). (29) 

The agreement between Y^o obtained using the sudden approximation, Eq. (pOj), and 
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the usual freeze-out calculation, Eq. (|), suggests that the sudden approximation is a 
reasonable point of departure for a first look at this phenomenon. Note, however, that 
the accuracy of the sudden approximation when self-annihilation is the principal photino 
equilibration mechanism does not guarantee that it is an equally good approximation 
when interconversion is the important process. The Boltzmann equation when photino 



self-annihilation dominates can be written 11 



dY _ -x{a~~\v\)s ^^^^ ,^2 
dx ~ Him) ^ ^Q^' 



{y'-yL\ (30) 



where Yeq{x) has the form given in Eq. ^ and H{m) = 1. 67 ^ /rripi. This is to be 
contrasted with the analogous expression when interconversion dominates: 

Here, is the equilibrium pion to entropy ratio: 

Y^{x) = 0.145(3/2)(2/10)(r^)x^/2exp(-r^x). (32) 

We have introduced = m-^/m = O.lTS/ig^, and included the factor (3/2) because the 
pion has three flavor x spin degrees of freedom in comparison to the photino's two. The 
difference in these forms, in particular the much weaker exponential dependence on x for 
compared to Yeq, is largely responsible for the shallower slope of the interconversion 
and inverse-decay curves as compared to the self-annihilation curves in Figs. 1 and 2. 
This shallower slope means that the quality of the sudden approximation in this case 
is inferior to the self-annihilation case, but probably not significantly in comparison to 
the large uncertainty due to our present poor knowledge of the cross sections. Closer 



examination of this question is in progress [|T^] 



Now we proceed using the sudden approximation. Given x*, we wish to determine 
fl-h"^. It is, of course, a very sensitive function of x^^: VL~h^ = 2.25 x 10^ [/ig] Y^o = 
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Fig. 3: ^~h^ as a function of assuming the photino stays in equilibrium until and immediately 
decouples (the sudden approximation). 

Tabic III: The value of fl~h^ assuming sudden freeze out aX x = x^,. fl~h'^ = 1 occurs around x^ ~ 20, 
and ^-h? = 10^^ around x* = 25. In the Table we have taken [i^ — 1. 
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47 


7.6 
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0.5 


0.2 


0.03 


0.01 


0.004 


0.0007 


10-6 



6.5 X 10^ [/ig] xj exp(— X*). The dependence of VL~Y? upon is shown in graphical 
form in Fig. 3, with specific values presented in Table 111. 

Since the age of the universe restricts fi-Zi^ to be smaller than unity, x* must be 
larger than 20. In order for the relic photinos to be dynamically interesting in structure 
evolution ^~h? must be larger than 10^^, which obtains for x^, < 25. Photinos would 
dominate the mass of the Universe if > 0.03,|^ which would result if x,, = 24. For 
Q~ = 1 and h ~ 1/2, x=k must be about 22. Thus we can summarize interesting values of 
^■^Nucleosynthesis bounds the contribution from baryons to be about flsh^ < 0.03. 
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1.2 1.4 1.6 1.8 2 2.2 

r=M/m (gluino-to-photino mass ratio) 

Fig. 4: Assuming 7 freeze out is determined by ^-BP conversion, the figure shows as a function of r 
the values of [mIms^C*] required to give the indicated values of Qr-h?. The uncertainty band is generated 
by allowing /xg to vary independently over the range 0.5 < /ig < 2. 



X* by 



< 25 

a;* <24 
20 < < 23 

X.. < 20 



^^^^ ~ 10~^; dynamically interesting role for photinos 

^7fh^ ^ 0.03; photinos dominate baryons 

fl-h"^ ~ 0.9, photinos are the dark matter and flTOT=i 



(33) 



^^^^ ~ 1; disallowed by age arguments. 
Now in turn, is exponentially sensitive to r = M/m, so limits to the contribution 
to the present density from 7 will be a sensitive probe of r. 

Prom Figs. 2 and 3, we see that for the canonical choices fi^ = /Js = A = B = C = 1, 
either the interconversion process, or decay-inverse decay is the last photino reaction 
to be of importance. It is impossible to say which one because of the uncertainties in 
the computation of the cross section and the decay width, so we shall consider both 
possibilities in turn. 

If interconversion determines the relic abundance and we make the sudden approxi- 
mation then we can determine as a function of the unknown parameters. Such a 
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r=M/m (gluino-to-photino mass ratio) 



Fig. 5: Assuming 7 freeze out is determined by decay /inverse decay, the figure shows as a function 
of r the values of [/x|/x^^-B] required to give the indicated values of 

graph is given in Fig. 4. Prom the graph we see that Q.-h?' < 1 can result for r — 2.2 if 
we allow fxlfXg'^C to be as large as 10^. We also see that a dynamically interesting value 
of can result for r as small as 1.2 if //|/x^^C = 10~^, although if the interconversion 
rate is suppressed this much, it is likely inverse decay would govern freeze out. 

A similar calculation can be made assuming that inverse decay is the last operative 
reaction depleting the photinos. The result of such an analysis is shown in Fig. 5. For 
r > 1.4 the behavior of the curves are similar to those in Fig. 5, but for small r the effect 
of phase-space suppression becomes important. 

In either case, the conclusion is that for r as large as 2.2, it is possible to have 
il-h'^ < 1; with our "central" choice of parameters, [j^lfXg^B] = [/ig^^/i^^'C] = 1, r must 
be less than 1.8 in order for fl-h"^ ^1- A value of r as small as 1.2 may result in 
fl-h"^ > 10~^; again with the central choice of parameters the limit is r > 1.6. 

Although it apparently is not important for realistic parameters, we mention a pos- 
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sible special role for the S*", udsg, the lightest baryon containing a gluino. Since the 
has a non-zero baryon number, its abundance is not given by Eq. (|^) at low temperature 
because of the non-zero baryon number of the Universe. So long as the strong inter- 
actions are maintaining equilibrium between nucleons and S'^'s, its abundance should 
be nso ~ nivexp[— (M50 — miv)/T], where un is the nucleon abundance and mjv is 
the nucleon mass. Thus at very low temperature its abundance will be larger than the 
abundance, so the co-annihilation and interconversion processes 75*" —>■ KN and 
7A^ —>■ KS^ are a potential sink of 7s which in principle could help keep the 7 in equi- 
librium. However for realistic cross sections, this does not seem to be important at the 
relevant temperatures. Likewise, although at low enough temperatures there are more 
nucleons than pions so that T~j^_^j^oj^ is larger than r~^_^^o^, freeze out has already 
occurred before the number density of nucleons begins to dominate that of pions. 



VI. SUMMARY AND CONCLUSIONS 



We have studied the reactions important for the decoupling and freeze-out of photinos 
having mass m less than about 1.5 GeV. We have found that it is crucial to include 
the interactions of the photino with the R^, the gluon-gluino bound state whose mass 
M is expected to lie in the range 1 to 2 GeV. The has strong interactions and 
thus annihilates extremely efficiently and stays in thermal equilibrium to much lower 
temperatures. In this circumstance, photino freeze-out occurs when the rate of reactions 
converting photinos to -R^'s falls below the expansion rate of the Universe. The rate 
of the 7-/2° interconversion interactions which keep photinos in thermal equilibrium, 
(77r < — > R^tt) or decay/inverse decay (771 < — > depends on the densities of 
photinos and pions, rather than on the square of the photino density, as is the case 
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for the self-annihilation process. For photinos of the relevant mass range (m ~ 800 
MeV), the pion abundance is enormous compared to the photino abundance. Therefore 
the photinos stay in equilibrium to much higher values oi x = m/T than they would 
if self-annihilation were the only operative process, resulting in a smaller relic density 
for a given photino mass and cross section. We find using the sudden approximation 
that light photinos are cosmologically acceptable for a range of 1.2 < r = M/m < 2.2. 
Within this range, if 1.6 < r < 2, the photinos are an excellent dark matter candidate. 
The precise range of r for which the photino accounts for the cold dark matter may 
shift when the sudden approximation is improved and cross sections are better known. 
However the general conclusion is robust: light photinos can account for the dark matter 
of the Universe for a suitable value of r, which is consistent with theoretical predictions 
in an attractive class of suSY-breaking mechanisms 0. 

Since 7-i?° interconversion governs freeze-out, the usual relation between VlK^ and 
the relic's annihilation cross section [1^ is not valid. If inverse decay is the operative 
process, then there is no direct prediction for the 7 scattering cross section on matter.[^ If 
77r < — > R^TT is the operative process, a quantitative solution of the Boltzmann equations 
can be used to infer its cross section. It will be significantly smaller, more-or-less by a 
factor n~{x^) /n.,^{x^)^ than the conventional cross-section used in planning relic detection 
experiments. 

Direct detection of low-mass relic photinos is more difficult than detection of high- 
mass (say m ~ 50 GeV) photinos. In addition to the low cross section mentioned above, 
the average energy deposition is {E) = m?MT{v'^) /{m + MtY where Mt is the target 

Since the short-distance dynamics entering the matrix element for —> is the same as for the 
scattering reaction — > R'^N, these could in principle be related. At this time however we do not have 
sufficient control of the hadron physics involved to make a quantitatively accurate theoretical prediction 
of the cross sections from the R^ lifetime. 
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mass. Thus existing and planned experiments using relatively heavy targets are not well 
adapted to this search. On the positive side, our photino is more likely to have spin- 
independent couplings to nucleons than expected in the conventional picture [0. This 



is because in the suSYbreaking mechanism which leads to the light photino and gluino 
under discussion here, the off-diagonal terms in the squark mass-squared matrix can be 
comparable to the diagonal terms.[3 



Indirect detection via annihilation of gravitationally concentrated photinos |l^ , [T9[| , for 
instance trapped in the Sun, is unlikely. Because they are low-mass WiMPS, evaporation 
is much more efficient than in the high-mass case, and they do not concentrate sufficiently. 
(And, of course, the cross section is smaller than conventionally supposed.) 

We also note that if the is stable, there will be a relic abundance of them, with 
an abundance relative to baryons determined by Mgo and xs = M50/T5, where T5 is 
the temperature of freeze-out. The 5° mass is expected to be 1.5 — 2 GeV, so let us 
define M50 = l.SyUi.s GeV. Then 

M50 — UIn' 



ns 1 /M5o\3/2 
— = 7 exp 



1 ^1.5/ii.5GeV\^/2 
4 



T 

iil^j exp[-(l-0.6/^,,)x,], (34) 

where the factor 1/4 accounts for the fact that the 5*° is a spin-0 state and comes in just 

one flavor, whereas there are 4 spin x flavor degrees of freedom for the baryons. The 

self-annihilation cross section should be comparable to that of the R^, so ignoring the 

difference between and masses, xs ~ rx^ji, where XfiR is the value of x at which 

r{R^R'^ X)/H = 1. From Figs. 1 and 2 we see that rxnR ~ 45, giving us/ub ~ 

7 X 10~^ for /X1.5 = 1, and smaller for larger /ii.5. Since the S'°'s are strongly interacting, 

even this small an abundance may be detectable. They will be more gravitationally 

^^See Ref. Q for aUowed ranges of the parameters determining the squark mass-squared matrix, 
H, tan/3, and Mq. 
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concentrated than standard wimps of comparable mass because they dissipate energy 
through their strong interactions, ahhough they do not form atoms or bind to nuclei]^. 

What, then, is the strategy for testing the proposal that photinos with mass less than 
or about 1 GeV constitute the cold dark matter of the Universe? Of course if an in the 
1 to 2 GeV range could be excluded by laboratory searches, our suggestion for the dark 
matter would be also excluded. Assuming though that these particles are discovered, 
knowledge of experimentally accessible properties of the photino and (in particular, 
their masses, the R^ lifetime, and the cross section for R^N —>■ jN) coupled with detailed 
numerical analysis of the freeze-out process, will allow a much more accurate prediction of 
the relic abundance than has been possible here. Since the relic density is exponentially 
dependent on r, which will one day be well measured, an accurate quantitative test of 
this idea will eventually be possible. 

In the meantime, theoretical work can elucidate the viability of this proposal. In the 

SUSY-breaking mechanism relevant to this scenario the parameters fi, tan/5, and Mq, 

which determine the photino and gluino masses, are highly constrained [Q. Relatively 

soon even more accurate predictions for the photino and gluino masses will be possible, 

narrowing the possible range of photino masses corresponding to any allowed gluino 

mass. With use of lattice gauge theory, it should be possible to compute the mass 

range corresponding to a given gluino mass, and thus to determine a spectrum of possible 

r values. Lattice gauge calculations can also in principle determine the matrix element 

for i?*^ — i> 77r, given the squark masses (which are fixed by the same unknown parameters 

determining m~ and m~, at least with minimal SUSY breaking), and determine the masses 

of the other i?-hadrons, which would help in estimating the cross sections. For instance 

knowledge of the mass of the Rt^ would allow one to better model o"/jo^. With more 

^^If they were stable and could bind to nuclei, they would have been detected in rare isotope searches||l| , 
so that possibility is excluded. 

29 



accurately fixed inputs, a full numerical solution of the coupled Boltzmann equations 
would be justified. 

Therefore, the most important next steps are: 

1. Look hard for i?-hadrons and other new particles predicted by this scenario. Planned 
kaon experiments may be able to establish evidence for the and possibly mea- 
sure the its lifetime and mass, as well as the mass of the photino f^. 

2. Do a better job fixing the parameters of the underlying theory, as well as calculating 
the photino mass produced through radiative corrections. 

3. Use lattice gauge theory to calculate the mass and check other predictions of 
this scenario such as the origin of the t]' mass 



4. A more complete treatment modeling the photino freeze-out is necessary |T^. An 
immediate question to address is the quahty of the sudden approximation used here. 
When interconversion is the dominant process, the equation governing the evolution 
of the photino density has a somewhat different form than in the self-annihilation 
case, for which the quality of the sudden approximation is well-established. 

5. Obtain detailed predictions for the low energy 7- nucleus cross sections expected in 
this scenario, and find effective detection techniques for light photinos. 

At the very least we have shown that until the value of r is demonstrated to be larger 
than about 2.2, light photinos are cosmologically acceptable. At best, we have described 
the scenario for the production and survival of the dark matter of the Universe. 

While there is no shortage of candidates for relic dark matter particle species, this 
proposal extends the idea that photinos may be the dark matter to a previously excluded 
mass range by incorporating new reactions that determine the photino relic abundance. 
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If this scenario is correct, direct and indirect detection of dark matter might be even 
more difiicuh than anticipated. However the scenario requires the existence of low-mass 
hadrons, which can be produced and detected at accelerators of moderate energy. Thus 
particle physics experiments will either disprove this scenario, or make light photinos the 
leading candidate for dark matter. 
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